Direct observation of quantum phonon fluctuations in a one dimensional Bose gas 
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We report the first direct observation of collective quantum fluctuations in a continuous field. Shot- 
to-shot atom number fluctuations in small sub-volumes of a weakly interacting ultracold atomic ID 
cloud are studied using in situ absorption imaging and statistical analysis of the density profiles. 
In the cloud centers, well in the quantum, quasicondensate regime, the ratio of chemical potential 
to thermal energy is ^/ksT ~ 4, and, owing to high resolution, up to 20% of the microscopically 
observed fiuctuations are quantum phonons. Within a non-local analysis at variable observation 
length, we observe a clear deviation from a classical field prediction, which reveals the emergence of 
dominant quantum fluctuations at short length scales, as the thermodynamic limit breaks down. 

PACS numbers: 



At temperature T = 0, classical thermodynamics pre- 
dicts the complete absence of excitations, however, due 
to the Heisenberg uncertainty principle, quantum ob- 
servables never fully come to rest. Quantum, i.e. vac- 
uum fluctuations are key to the understanding of quan- 
tum electrodynamics effects as fundamental as sponta- 
neous emission [T], the Lamb shift or Casimir-Polder 
forces near surfaces [3H5], but also. Hawking radiation 
near black holes [6], quantum phase transitions jTj, etc. 
In matter fields, quantum fluctuations govern the cor- 
relations properties at low temperature. They cause 
quantum depletion in Bose-Einstein condensates, bring- 
ing corrections to their equation of state (EoS) [5111], and 
are dramatically enhanced in reduced dimensions [10' . In 
ID systems, they destroy long range order and prevent 
Bose-Einstein condensation even at T = |11] . 

So far, quantum fluctuations in continuous flelds have 
been detected only indirectly, from their macroscopic 
consequences at the thermodynamic scale [SJ H], while 
their microscopic observation has remained elusive |12i 
I13| . Only recently, density imaging of ultracold atomic 
clouds has allowed to microscopically image a quantum 
fluctuating fleld (more precisely, its modulus) in discrete 
lattice systems |!l4j. In this Letter, we report microscopic 
observation of quantum fluctuations for the flrst time in 
a continuous fleld. 

Statistical analysis of in-situ density fluctuations has 
indeed become a prominent tool of investigation. In di- 
mension D, the variance of atom number {SN'^) in a suf- 
flciently large volume is the same as for an inflnite 
system in thermodynamic equilibrium at the local den- 
sity n and chemical potential fi |15l I16| : 

{6N^) - A''kBT{dn/dfi)T, (1) 

where fc^ is the Boltzmann constant, and {dn/d^)T /n^ 
the isothermal compressibility, derived from the EoS 
n{p,T). This thermodynamic regime has allowed to ob- 
serve bosonic bunching |15} I17| and fermionic antibunch- 
ing jlHlIT^ in ideal gases. In weakly repulsive Bose gases, 
the suppression of density fluctuations, which deflnes the 



quasicondensate regime, was detected directly in ID |17] 
and 2D [20 . Recent thermodynamic studies addressed 
the phase diagram for quasicondensation in ID and elon- 
gated 3D Bose gases [21] , the universality of the 2D Bose 
gas [20,, and, in optical lattices, the superfluid to Mott 
insulator transition [22]. In ID Bose gases at fc^T < /i, 
subpoissonian density fluctuations, i.e., antibunching [25] 
provided thermodynamic evidence of the quantum quasi- 
condensate regime [241 125| . 

Thermodynamic measurements along Eq. [l] require 
that the relevant excitations under observation have oc- 
cupation numbers ^ 1, so that thermal fluctuations 
dominate, and fluctuations are classical, i.e. proportional 
to T, as in Eq. ^ with no quantum contribution. This is 
valid only in the thermodynamic limit 

Lohs > la, (2) 

where Z J is the thermal correlation length of density fluc- 
tuations, and Lobs is the length scale at which the system 
is probed. If Lobs < , measured fluctuations can devi- 
ate from Eq|l] |16l I2n] . and, if < 1, quantum fluctua- 
tions dominate. In this Letter, we analyze density fluc- 
tuations in ID Bose gases with /i/fc^T ~ 4, reaching the 
crossover regime ~ Lobs- Here, quantum fluctuations 
have a sizeable contribution. Still, quantum and thermal 
non-thermodynamic effects cancel almost exactly, and 
Eq. [l] apparently holds. However, quantum fluctuations 
are non-extensive |16|, I23| . scaling only logarithmically 
with Lobs, while thermal fluctuations scale linearly (as 
in Eq. [l] with Lobs — A). Hence, truly thermodynamic, 
i.e., purely thermal fluctuations, are observed only for 
Lobs ^ l]: ■ By varying Lobs, we monitor the breakdown 
of the thermodynamic limit at small Lobs as a clear devi- 
ation from a classical fleld theory that ignores quantum 
fluctuations. In the line of [13J, this criterium proves our 
microscopic observation of quantum fluctuations. 

Our experiment uses ®^Rb atoms in micro-magnetic 
traps on an atom chip. The transverse and longitu- 
dinal trapping frequencies are luj_/2tt ~ 3.3 kHz and 
a;||/27r = 5.5 Hz respectively. After forced rf evapora- 
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FIG. 1: (Color online). Microscopic fluctuations and ther- 
modynamic theories, at T = 18 nK (b) and T = 4.7 nK (a, 
c). (a) : Typical absorption picture (optical density), (b, c) 
Normalized fluctuations and average profiles (insets). Theo- 
retical curves : classical shot noise (dotted), predictions from 
Eq.[l]and the ideal Bose gas (short-dashed), quasi-condensate 
(long-dashed), and Yang- Yang EoS (solid). Grey areas indi- 
cate in the profiles indicate the quantum quasi-condensate 
regime, the dark-grey area in the T = 4.7 nK profile is the 
zone f2 where the analysis of Fig. 4 is carried. 
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FIG. 2: (Color online). Quantum vs thermal fiuctuations in 
ID quasicondensates. (a, b) Fluctuation regimes for thermal 
(a) and quantum (b) quasicondensates. (c) (solid) and Sl_ 
(dashed) in Bogoliubov theory for various ksT / fi. (d) ^'^^^(z) 
for 7 = 0.01. 



live cooling and thernialization for 800 ms, an absorption 
picture is recorded on a CCD camera. After hundreds of 
realizations, fluctuations in the density profiles are an- 
alyzed as detailed in p3]. For each profile and pixel of 
length A = 4.5 /im, we extract the atom number fiuc- 
tuation 6N ~ N — (N) where (N) — nA is the mean 
atom number. True atom number variances {SN'^) are 
inferred from measured variances {5N'^)m using the ther- 
modynamic relation {6N'^)m = K2{5N'^), where K2 is a 
reduction factor due to the finite rms imaging resolution 
6. Assuming a gaussian imaging response, S is deter- 
mined precisely from the correlations measured between 
neighboring pixels |15l [27] . 

To introduce our data, we show in Fig. 1 a typical pic- 
ture (a), average profiles, and relative fluctuations, for 
clouds deep in the ID regime, with ksT/hoj^ =0.11 (b) 
and 0.03 (c). As in [151 1^ . T is measured accurately by 
fitting the fiuctuations in the cloud centers with Eq. [l] 
and the Gross-Pitaevski (G-P) quasicondensate EoS in 
the 1D-3D crossover /i = fujj^{y/l 4- 4na— 1) [2H] (dashed 
curve, see also Fig. 3), where a = 5.3 nm is the 3D scat- 
tering length. At low densities, fiuctuations are super- 
poissonian, i.e., exceed the shot noise (dotted line). They 
reach a maximum at the transition to the quasiconden- 
sate regime, around the density riqc — [m(fcBT)^)//i^(7]^/'^ 
|21| . in good agreement with the exact Yang- Yang solu- 
tion for the ID Bose gas |29j (solid curve) [31]. At higher 
densities, for fi > ksT (light grey areas in profiles), 
quantum fiuctuations start dominating thermodynami- 
cally and fiuctuations turn subpoissonian, as explained in 
|23j . In this Letter, the important region is the cloud cen- 
ter at T = 4.7 nK (zone n), where {SN^)/{N) ~ keT/fi 
reaches 0.26, and where we detect quantum phonons mi- 



croscopically, as explained below. 

To analyze fiuctuations in the cloud centers, we use 
Bogoliubov theory, valid for weakly interacting quasi- 
condensates f5D], and the ID G-P EoS fi = gn, where 
g — 2fiLd±a is the ID coupling constant. This is appropri- 
ate since the dimensional correction to the G-P prediction 
is only 10% in zone 57, and the ID interaction parame- 
ter is 7 = gm/h^n < 0.05. Bogoliubov excitations have 
energies e/c — h'^k'^ /2rn^/Y^^4jW^, where ^ = h/y/nifi 
is the healing length, and thermal occupation numbers 
Uk = l/(e"''/'=«'^ - 1). Noting fk = l/^/m/k^, the 
spectrum of density fiuctuations is 

{Snl)=nSk = niS^ + Sl), (3) 

where — fk and 5j = 2/fcni. are the quantum and 
thermal static structure factors, both plotted in Fig. 2c. 
Note that, in Sk = 2fk(nk + 1/2), S'^ is the exact analog 
of the zero-point energy term in the harmonic oscillator 
spectrum. We also show in Fig. 2d the second order 
correlation function g^'^\z) = 1 + / ||e*'^^(S'fc — 1). 

Three competing terms determine the fiuctuation 
regimes in quasicondensates. For fc^ ^ 1, — 1 is 
the trivial autocorrelation "shot noise" |15[ [T7] of un- 
correlated particles, for which g^'^\z) = 1. At smaller 
fc, repulsive interactions give a negative contribution to 
5fc - 1 and to g^'^\z) - 1 [23]. Finally, in Sl , statistical 
bosonic bunching brings a positive contribution. Let us 
now focus on the quantum and thermal density correla- 
tion lengths ZJ? and , respectively defined as the inverse 
widths of s2 and S'J. For quantum fiuctuations, equally 
present in all modes fc, one has = £, [ST. For thermal 
fiuctuations, two cases need be distinguished. 

In thermal quasi- condensates {k^T 3> ^, see Fig. 2a), 
5 J decays fast for kl^ ^ 1, so one has a single Ic = 
= [251 [3T] . However thermal fiuctuations domi- 
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FIG. 3: (Color online). Quantum vs thermal contributions in 
the measured fluctuations {5N'^)m (same data as in Fig. 1). 
Dot-dashed lines are the spectral ID Bogoliubov predictions 
{5N'^)'i (low-lying green), {5N'^)l, (middle red) and {5iV^)|^' 
(upper blue). Thermodynamic predictions as in Fig. 1 are 
here rescaled by K2 = 0.34 (a) and 0.41 (b) for 5 — 3.5 and 
2.8 fim, respectively. The square in (b) is the zone Q,. 



nate, i.e. rik ^ 1, for length scales L ^ 1/k ^ AdB, 
where AdB = y/^jrWJmkBT is the thermal de Broglie 
wavelength. Three regimes are thus present |16| . For 
L ^ Ic, thermal phonons dominate : this is the thermo- 
dynamic regime, with superpoissonian fluctuations. For 
'^dB L <^ Ic, thermal particles dominate, but correla- 
tions are partly lost [37]. Finally, for L <^ AdB, only the 
poissonian quantum shot noise is detected. 

In quantum quasicondensates (ksT ^, see Fig. 2b), 
noting c — y/jjjm the speed of sound, Uk and S"^ are 
both suppressed at length scales lower than [5H] 



he 



= i9 



(4) 



which, to our knowledge, has been proposed only recently 
|13LfTBl[5?| . The familiar regimes of poissonian shot noise 
and thermodynamic thermal phonon fluctuations (now 
subpoissonian) are still present. However, in the range 
<C i ^ ZjT, one now has a crucial new regime (high- 
lighted in Fig. 2b), where <C 1 and quantum phonons 
govern the physics. 

To compare Bogoliubov theory to our data, we use 
the imaging transfer function Uk ~ ^ — 



f) ^i-cos(fcA)) „-k^S^ 



|23| . obtained from our gaussian optical response model 
[15t[27]. and we compute, for j — Q,T, 



{SN 



2\i 



dk 
2^ 



nSlUk- 



(5) 



Uk is peaked at fc = 0, and its inverse width is Lobs ^ 
max{(5, A}. The thermodynamic limit Eq. [2] is thus 
equivalent to measuring {SN'^)„i = K2{N)So [ini E7] . 
i.e., to only probe the contribution 5*0 = Sq of ther- 
mal phonons, always proportional to T (see Fig. 2c). In 
other words, since Sq = 1 + f dz{g'-^\z) — l), a thermody- 
namic observation probes only the integral of g^'^\z) — 1, 
without resolving its microscopic details (see Fig. 2d). 

In Fig. 3, we compare the measured fluctuations 
{SN'^)m to the spectral ID Bogoliubov predictions 
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FIG. 4: (Color online). Fluctuations at variable length scales 
in the cloud center (zone Q) at T = 4.7 nK. (a) Typical profile 
TV (triangles), reference profile (A*') (circles) and fluctuations 
SN (squares), for three pixel sizes Lpix. (d) (SN^)/ (N) in zone 
fl compared to full Bogohubov (FBM) prediction {SN'^f^^ /k2 
(solid), classical field (CFM) prediction (5N'^)'^/k2 (dashed) 
and thermodynamic G-P prediction (dot-dashed) at Tjq = 
5.0 uK. Inset : same data fitted with FBM at T^o = 5.0 nK 
(sohd) or CFM at T^^ = 5.5 nK (dashed). 



{5N^)g and {SN^)l„ and their sum {SN^)2 (dot-dashed 
lines). At low densities, one has l'^ ^ Lobs, which ex- 
plains that ((5iV^)^ follows the shot linear noise predic- 
tion (see inset to Fig. 3a). At higher densities, {SN^)^ 
is reduced by repulsive interactions [33] • At T = 18 nK, 
the ratio {SN'^)'^/ {SN"^)^ is 7% and cannot be resolved, 
as in [25. However, at T = 4.7 nK, it reaches 20% in the 
cloud center, which now exceeds the experimental uncer- 
tainty. Thus, the contribution of quantum fluctuations 
is here sizeable, i.e. non-negligible, in each pixel, in con- 
trast to the thermodynamic regime. Yet, {6NY^ is still 
in good agreement with the ID thermodynamic Yang- 
Yang prediction (solid line) [30] • This is because and 
5j have opposite slopes ±fc^/2 at small k^ for all T (see 
Fig. 2.c), so that non-thermodynamic (i.e., finite k |16| ) 
quantum and thermal contributions cancel each other at 
first order. 

To confirm our detection of quantum fiuctuations, we 
turn to a non-local analysis. Since quantum fiuctuations 
scale only logarithmically with Lobs [16, 23_, thermal fiuc- 
tuations always dominate for Lobs ^ and fiuctuations 
are well described by a classical field model (CFM) that 
ignores the quantum term S*^ [13J, On the other hand, a 
CFM is expected to fail for Lobs ■ To check this, we 
focus on the zone il, i.e., the 3 central bins 39 < {N) < 54 
at r = 4.7 nK, containing 58% of the atoms (see inset 
to Fig. Id and Fig. 3b). There, the profile is the fiattest, 
many data points are available, and error bars are the 
smallest. The cloud center is indeed the most reliable 
fraction of the data. 

To vary Lobs, we merge the imaging pixels in macro- 
pixels of variable size npix = 1 to 8, and compute fiuctua- 
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tions accordingly (see Fig. 4a and |17]). Figure 4b shows 
{SN^)/{N) obtained in the zone fl, for Lpix = Anpix = 
4.5 to 36 /im. The key criteriuni. as in [13_, is to com- 
pare the data either to a full Bogoliubov model (FBM) 
that includes the quantum term (solid line) or to a CFM 
that ignores it (dashed line). At large ipix (— ^obs), both 
models converge as expected towards the thermodynamic 
ID G-P prediction (dot-dashed line), with a logarithmic 
vanishing of the quantum contribution. At short Lpix, 
the data clearly deviate from the CFM prediction, which 
displays a noticeable dip. Note that the gap saturates at 
small Lpix, due to the finite optical resolution that cuts 
off high k fluctuations, which conflrms that shot noise is 
irrelevant in our data. 

The theoretical predictions in Fig. 4b are computed 
at T*j^ = 5.0 nK, that we obtain fitting the FBM to the 
data, with T as only free parameter (see inset to Fig. 4b). 
This fit has an rms deviation of only 1.0%, i.e., much 
less than the plotted error bars [H] and is thus a very 
accurate ID thermometry |42]. On the contrary, the fit to 
the CFM, yielding T-fjj = 5.5 nK, has a strong systematic 
error and an rms deviation of 5% (see inset to Fig. 4b). 
This clear breakdown of the CFM proves our observation 
of quantum fiuctuations, and reveals the emergence of 
dominant quantum phonons at short distances. 

As for length scales, in the zone fi, l'^ — 0.6 /im and 
= 2.2 nm, while Lobs is determined by A = 4.5 /im 
and 6 = 2.8 /im. Thus, Lobs ~ Ic ^-^d this explains 
our unprecedented observations. On the contrary, in 
|12| . at r = 33 nK, one has l]^ — 0.63 /im, and 
Lobs > A = 10 iim, i.e.. Lobs ^ ^if • All observations in 
|12j were thus well in the thermal fiuctuations regime, in- 
distinguishable from a CFM prediction [T3] . In the quan- 
tum regime z ^ , the decay of the first order (phase) 
correlation function (/'^-"(z) is algebraic, whereas, in the 
thermal regime z ^ l]^ , it is exponential, over a thermal 
phase correlation length lT = I'^/^/j » [TTl [551151]. 
The misinterpretation in [121 came from identifying the 
quantum regime with z l'^. 

In summary, we have reported the first microscopic 
observation of vacuum fiuctuations in a continuous field, 
using a non-local analysis that reveals a clear deviation 
from a classical field theory. Our observation of emerging 
dominant quantum phonons is a first microscopic insight 
into the regime of quasi-long range order, i.e, algebraic 
decay of g(-'^)(z), in the ID Bose gas |TT]. We also demon- 
strated the possibility of imaging vacuum phonon fiuctu- 
ations in single density pictures, like Fig. la. By further 
reducing T and Lobs, one could monitor the full crossover 
from thermal to quantum fiuctuations [27]. Dark soli- 
tons, which are defects localized over a length scale ^ l^ , 
could also be used as sensitive probes for the microscopic 
physics of quantum fiuctuations jSH |35] . 
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